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Models with dark energy decaying into dark matter have been proposed in Cosmology to solve 
the coincidence problem. We study the effect of such coupling on the cosmic microwave background 
temperature anisotropies. The interaction changes the rate of evolution of the metric potentials, the 
growth rate of matter density perturbations and modifies the Integrated Sachs- Wolfe component of 
cosmic microwave background temperature anisotropies, enhancing the effect. Cross-correlation of 
galaxy catalogs with CMB maps provides a model independent test to constrain the interaction. 
We particularize our analysis for a specific interacting model and show that galaxy catalogs with 
median redshifts z m = 0.1 — 0.9 can rule out models with an interaction parameter strength of 
c 2 ~ 0.1 better than 99.95% confidence level. Values of c 2 < 0.01 are compatible with the data and 
may account for the possible discrepancy between the fraction of dark energy derived from WMAP 
3yr data and the fraction obtained from the ISW effect. Measuring the fraction of dark energy by 
these two methods could provide evidence of an interaction. 



INTRODUCTION 



Measurements of luminosity distances using supernovae type la (SNIa) of the cosmic microwave background 
(CMB) temperature anisotropies with the WMAP satellite Q, large scale structure 0,0], the integrated Sachs- 
Wolfe effect [5|, @ , and weak lensing Q , indicate that the Universe has entered a phase of accelerated expansion 
-see Q for recent reviews. This acceleration is explained in terms of an unknown and nearly unclustered 
matter component of negative pressure, dubbed "dark energy" (energy density, p x ), that currently contributes 
with about 75% of the total density. The remaining 25% is shared between cold dark matter (p c ~ 21%), 
baryons (pt ~ 4%), and negligible amounts of relativistic particles (photons and neutrinos). The pressureless 
non-relativistic matter component redshifts faster with expansion than the dark energy, giving rise to the so- 
called 'coincidence problem' that seriously affects many models of late acceleration, particularly the ACDM 
model: "Why are matter and dark energy densities of the same order precisely today?" Q. To address this 
problem within general relativity one must either accept an evolving dark energy field or adopt an incredibly 
tiny cosmological constant and admit that the "coincidence" is just that, a coincidence that might be alleviated 
by turning to anthropic ideas 10]. (It should be mentioned, however, the existence of proposals in which a 
vacuum energy density of about the right order of magnitude arises from the Casimir effect at cosmic scales 
-see [ll[ and references therein). 



One way to address the coincidence problem, within general relativity, is to assume an interaction (coupling) 
between the dark energy component and cold dark matter such that the ratio r = p c j p x evolves from a 
constant but unstable value at early times (in the radiation and matter dominated epochs) to a lower, constant 
and stable attractor at late times well in the accelerated expansion era [12|, Il3| . Before the late accelerated 
expansion was observed, Wetterich introduced interacting quintessence models to reduce the theoretical high 
value of the cosmological constant [l4j]. Later, these kind of models were rediscovered, sometimes (but not 
always) in connection with the coincidence problem -see, e.g. [HI]. Other solutions (known as f(R) models) 
require the Einstein-Hilbert action to be modified fig . 



Interacting quintessence models are testable since they predict differences in the expansion rates of the 
Universe, the growth of matter density perturbations, the Cosmic Microwave Background (CMB) temperature 
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anisotropics and in other observables. In this paper we shall demonstrate that, as the Integrated Sachs- Wolfe 
effect measures directly the growth rate of matter density perturbations, can be used to detect variations on 
the evolution of the large scale structure with respect to the prediction of the concordance ACDM model. As 
a toy model, we shall consider a spatially flat Friedmann-Robertson-Walker universe filled with radiation, 
baryons, dark matter and dark energy such that the last two components interact with each other through a 
coupling term, Q = 3TLc 2 (p c + p x ). Thus, the energy balance equations for dark matter and dark energy take 
the form 

p c + 3Hp c = 3Hc 2 (p c + p x ) , p x + 3H{1 + w x )p x = ~3Hc 2 (p c + p x ) , (1) 

where w x = p x /p x < —1/3 is the equation of state parameter of the dark energy component and c 2 is a constant 
parameter that measures the strength of the interaction (it should not be confused with the speed of light 
that we set equal to unity). Derivatives are taken with respect to the conformal time, and 7i = a/a, with a 
the scale factor of the Friedmann-Robertson-Walker metric. To satisfy the severe constraints imposed by local 
gravity experiments [l7l Il8j , baryons and radiation couple to the other two energy components only through 
gravity. Our ansatz for Q guarantees that the ratio between energy densities, r, tends to a fixed, attractor, 
value at late times. It yields a constant but unstable ratio at early times. The details of the calculation can 
be found in [l3|, [2(J This result holds irrespective of whether the dark energy is a quintessence (i.e., 
— 1 < w x < —1/3), phantom (w x < — 1), fc-essence or tachyon scalar field [22| . It is worth to remark that the 
above ansatz implies an effective power law potential for the quintessence dark energy field at early times and 
an effective exponential potential at late times 

Beyond the physical motivation for the existence of a coupling between dark matter and dark energy, their 
interaction parameter c 2 needs to be constrained by observations. In [2(| we showed that the slope of the 
matter power spectrum at small scales decreases and the scale of matter radiation equality is shifted to larger 
scales with the interaction. CMB temperature anisotropies are equally affected if the interaction is present 
at decoupling. The luminosity distance of SNIa depends on the expansion rate of the Universe and can be 
used to test the existence of a DM/DE coupling. Comparison of model predictions with data showed that only 
WMAP 3yr data and the Sloan Digital Sky Survey (SDSS) data provided constraints on the value of c 2 [20, 21]. 
In particular, the SNIa data proved to be rather insensitive to the coupling. In this paper we show that the 
Integrated Sachs- Wolfe (ISW) component of the CMB temperature anisotropies is much more sensitive than 
the SNIa data at redhisfts z ~ 1. It measures the growth factor of linear matter density perturbations, that 
is strongly affected by the interaction and can effectively set constraints on Interacting Quintessence Models 
(IQM). Other experiments that also constrain the growth rate are discussed in (23|. In Section II we briefly 
summarize the formalism behind the measurements of the ISW effect by cross-correlating CMB data with galaxy 
catalogs. In Section III we describe the observational prospects and in Section IV we present our conclusions. 



II. THE INTEGRATED SACHS- WOLFE EFFECT 



When CMB photons in their way to the observer from the surface of last scattering, cross a time-varying 
potential well, experience the so called late Integrated Sachs- Wolfe effect. Photons are blueshifted (redshiftcd) 
when entering (leaving) a high dense region, and similarly for low dense regions. If the gravitational potential 
(f> evolves during a photon crossing like in models with a cosmological constant, then both effects will not 



compensate each other and the final energy of the photon will vary. The net effect will be [24 1 



« = _ 2 r-^, (2) 



To J dr 

The integral is taken with respect to conformal distance (or look-back time) from the observer at redshift z = 0, 
r(z) = L dz' /H(z') where the Hubble factor, H, obeys Friedmann's equation, 3H 2 = 8irG(p c + p x ). 



The evolution of the Newtonian potential will be computed numerically, but its qualitative behavior can be 
understood from the perturbed Einstein's equations. Once the Universe becomes matter dominated, we can 



3 



0.0 

-0.2 
-0.4 
-0.6 



1 .0 
1 00.0 





' ' 








(a) : 




\\ 


















^ \ 




- 




_ 




'\ \ 






\ \ 






\ \ 




- 


\ \ 


- 




\ \ 






\ \ 






\ \ 






\\ 












V 











0.8 
0.6 

0.4 



0. 



0.0 



1 .0 



0.1 



00 



1 



FIG. 1: Panel (a): Evolution of the background equation of state parameter, wb, as a function of redshift. Panel (b): 
Evolution of the Newtonian potential <j) for the wavenumber mode k = O.Ol/i Mpc -1 . Solid, dashed and dot-dashed lines 
correspond to c 2 = 0,0.01, and 0.1, respectively. In drawing the figures we have set w x — —0.9. 



write in the conformal Newtonian gauge [2fJ, \2l 



H(f>-3w B H 2 



Sifli 



(3) 



where the sum is over all matter fluids and scalar fields, Qi, Si are the density parameter and the energy density 
contrast, respectively, of fluid i. Derivatives are also taken with respect to conformal time. In non-interacting 
models, the gravitational potential changes when (a) the equation of the state of the background varies, and 
(b) the dark energy is smoothly distributed, i.e., it does not cluster on scales below the "sound horizon", set 
hy c 2 s x = (SP x /5p x ). At that point the potential will decay and generate an ISW effect (2(| [27|- During the 
matter dominated epoch, the matter sound speed is c 2 c = (SP c /Sp c ) = 0. If c 2 x = 1, then the dark energy 
fluid does not cluster on scales below the horizon [28| . whereby 5 X ~ 0. In this case, on subhorizon scales the 
gravitational potential evolves as 4> = Ci + C^a -5 / 2 , i.e., it essentially remains constant. However, in models 
with interaction, the CDM energy density does not scale as a~ 3 because of the coupling. During the matter 
dominated regime, the background equation of state parameter is wb = c 2 w x /(c 2 + 1) < and the growing 
mode in Eq. becomes cf> — C\ a rn with m < 3uib/4 < 0; that is, there are no growing modes and the 
potential <fi decays. In the period of accelerated expansion, in both interacting and non-interacting models, </> 
decays. In Fig. []Ji we plot the evolution of wb as a function of redshift. Solid, dashed and dot-dashed lines 
correspond to c 2 = 0,0.01, and 0.1, respectively, and we have set w x = —0.9. Notice that when c 2 ^ 0, the 
matter dominated period has ^ 0. In Fig. [Tp we depict the time evolution of the Newtonian potential (j> 
as a function of redshift for a mode of wavenumber k = 0.01/iMpc -1 . In the three cases, the amplitude of 
the potential is normalized to unity at horizon crossing. We solved the system of equations that describes the 
evolution of all energy components and the Newtonian potentials using a modified version of the CMBFAST 
code 29]. As the figure shows, when c 2 = 0, 4> remains constant during the matter dominated epoch, and it 
decays solely during the phase of accelerated expansion. By contrast, in interacting models the potential 4> 
decays even during matter domination. 



The ISW effect is not the only anisotropy generated by local structures. Secondary anisotropics could be 
generated due to non- linear evolution of the gravitational potential, usually associated with the gravitational 
collapse of small scale structures like clusters [3(| . But the relevant scales are those of clusters and superclusters, 
corresponding to angular scales of 5 — 10 arcmin, much smaller than those associated to the ISW effect. 
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FIG. 2: (a) growth rate of matter density perturbations, (b) growth rate relative to Einstein-de Sitter and (c) differential 
growth rate per unit of redshift for two cosmological models with w x — —0.9 and different dark energy fraction and 
interaction parameter. In all panels, solid lines correspond to £l x = 0.8, and dashed lines to £l x = 0.7. In each set and 
from top to bottom, the interaction parameter is c 2 = 0.1, 0.01, and 0. 



The ISW component cannot be measured directly, i.e., cannot be extracted from CMB data using frequency 
information since it the same frequency dependence that any other intrinsic CMB signal. Since at low 
redshifts is generated by local structures, the spatial pattern of ISW anisotropics will correlate with the matter 
distribution in the local Universe. The contribution to the radiation power spectrum due to the evolution of 
the gravitational potential at low redshifts is given by [3l|, [HJ 

Cl sw =< \a lm \ 2 >= - / k 2 dkP(k)lf(k) , (4) 
n Jo 

where 

fT2 (Tree JP 

I l (k)=m c , ^J dr 3l {kr) — . (5) 

In this expression, F = D + (z)(l+z) is the ratio of the growth of matter density perturbations in any cosmological 
model, D + (z), to that in a Einstein-de Sitter Universe, (1 + z)^ 1 ; Hq 1 is the present Hubble radius and ji{x) 
are spherical Bessel functions. In Figure [2] we plot the redshift evolution of D + , F and dF/dz. All the models 
have w x — —0.9 and solid lines correspond to Q x = 0.8 and dashed lines to fl x — 0.7. From top to bottom, each 
group of lines corresponds to c 2 = 0.1, 0.01, and 0, the latter pertaining to non-interacting models. The growth 
function was normalized to unity at z — 0, D+(0) = 1. 

Due to the ISW component, CMB temperature maps will correlate with tracers of the gravitational potential 
at low redshifts that could be constructed from galaxy catalogs. In particular, it will correlate with templates 
of the projected galaxy density along the line of sight. Assuming the distribution of galaxies follows that of 
the dark matter, the correlation of a template of projected galaxy density and a CMB map can be computed 
analytically (3l| . 



ISW-M 



^ J fc 2 dfcP(fc)// 5M/ (fc)/ ; M (fc), (6) 



with 



drW M {k ir )ji(kr), 



(7) 
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where W M (k,r) = b{k,r)nM{r) D + {r) is the window function of the tracer of the large scale gravita- 
tional potential, b(k, r) is the scale dependent bias existing between sources and the matter distribution as 
a function of radial distance r, tim(?") is the normalized rcdshift distribution of sources and D+(r) is the 
growth function of matter density perturbations. Similarly, the power spectrum of the galaxy catalog is given by 



One advantage of using Eq. © is that it separates the contribution coming from redshift evolution and from 
the spatial distribution of the large scale structure. Models of structure formation that give rise to the same 
matter power spectrum but differ in their growth rate, D + (z) (see Fig. [2b) , will also give rise to ISW effects of 
different amplitude. Consequently, the correlation with galaxy catalogs will also differ. If the interaction varies 
with time and is only significant at present, the shape of the matter power spectrum will be unaffected but 
the growth rate will not. Therefore, the ISW effect will be sensitive to the existence of the interaction through 
the influence of the latter on the growth rate of density perturbations. In the next section, we shall fix the 
matter power spectrum of all models with interaction to be the same as a model with c 2 = (no interaction) 
and the same cosmological parameters. We considered only power spectra that were compatible in amplitude 
and shape with the spectra derived from the SDSS and 2dGF galaxy catalogs [4J. Thus, differences in the 
correlation between templates of projected galaxy density and CMB data -Eq.©- will reflect differences in the 
growth ratio F only and will indicate how sensitive the ISW effect is to the redshift evolution of matter density 
perturbations. 



III. OBSERVATIONAL PROSPECTS 



After the release of WMAP lyr data, several groups looked for evidence of ISW effect by cross-correlating the 
WMAP data with templates built from different catalo gs: the NVSS radio survey [H, E3 , the HEAO-I X-ray 
Background map [5], the APM j^, the SDSS @, EE TH3], and the 2MASS galaxy surveys [H, El- These 
analyses rely on constructing templates from catalogs as tracers of the ISW signal. Cross-correlation of random 
realizations of CMB data and noise provides an estimate of the uncertainties. Catalogs differed in depth and 
sky coverage and different groups use different techniques: real space correlation, correlations in Fourier space, 
wavelets analysis, etc., and obtained positive detections with levels of significance of 2 — 4cr. These detections 
represent an independent evidence for the existence of dark energy, though, at the moment, the measurements 
are not accurate enough to discriminate between competing dark energy scenarios. 

In this section, we analyze the prospects that the ISW effect could be potentially sensitive to the growth rate 
of matter density perturbations and, therefore, to the dark energy-dark matter coupling. The effect of the 
interaction is to increase the amplitude of the cross-correlation of template and CMB data: the interaction 
makes the gravitational potentials decay faster which produces a larger ISW effect, making it easier to detect. 
If we trace the ISW effect with a template M that covers a fraction f s f. y of the sky, the cumulative signal to 
noise ratio is given by _32] 

g\2 'max rfjISW-M Y 

where C ; Ai is the power spectrum of the template M. The sum extends to the largest angular scale resolved by 
the template. The main source of confusion in ISW measurements comes from the intrinsic CMB anisotropies. 
In Eq. © we have neglected the instrumental noise of WMAP data and the shot noise of the template map. 

If the distribution of galaxies as a function of redshift is known, then the expected S/N ratio can be computed 
numerically. Following earlier practice, we model the redshift distribution of sources by the analytic expression, 
nu(z) — A{z/zq) 2 exp[— (z/zq) 3 ^ 2 ]. The constant zq measures the effective depth of the catalog and is simply 
related to the median redshift z m of the source distribution, z m = 1.412 zq [40H4l|. The normalization constant is 
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FIG. 3: Cumulative signal to noise ratio for a catalog with mean depth z m = 0.3. From bottom to top lines correspond 
to c 2 = (bottom thin-solid line), 10 -3 (dot dashed line), 0.01 (dashed line), and 0.1 (solid line). Template and data 
are assumed to have full sky coverage. 

fixed by setting J riM{z)dz = 1. In Fig.[3]we plot the cumulative S/N ratio for a cross-correlation of temperature 
data with a template constructed from a galaxy catalog with median redshift z m = 0.3, as a function of the 
multipole l m ax for 4 different interaction parameters. From bottom to top c 2 = (bottom thin-solid line), 
10 -3 (dot dashed line, almost coincident with the previous one), 0.01 (dashed line), and 0.1 (solid line). For 
simplicity we considered data and template to have full sky coverage, i.e., f s k y = 1. In Fig. [3] it is seen that 
increasing the strength of the interaction augments the S/N ratio reflecting the fact that, as the interaction 
generates larger ISW signal, it can be detected with higher statistical significance. 

Comparison between model predictions and observations can be done by means of the correlation function. The 
theoretical prediction is < T*M > {9) = £,(2/ + l)C{ sw ~ M P t (cos 0)/4tt. The cross-correlation of a template 
M constructed from a tracer of the large scale matter distribution like a galaxy catalog with a CMB map is simply 



where Mi is the projected galaxy density on pixel i and <Ji is the noise associated to each pixel on a CMB map. 
The measured correlation function in Eq. (fTOf is affected by masking, pixelization and the inhomogeneities 
within the catalog of the galaxy selection function across the sky. Further, the free electrons that reside in 
the deep potential wells of galaxy clusters contribute to the temperature anisotropics via the thermal Sunyaev- 
Zcldovich (SZ) effect (HJ, that dilute the ISW signature at zero lag Q. For example, templates constructed 
with the 2MASS catalog of galaxies are known to trace the thermal SZ component [43[ as well the ISW effect. 
To avoid a possible thermal SZ contamination, we will compare the mean between 4 and 10° of the predicted 
and measured correlation functions. In Figure 2] we plot the correlation function computed using the analytic 
galaxy redshift distribution given above. Theoretical estimates were done assuming the bias factor b to be 
constant. Dashed lines corresponds to models with no interaction (c 2 = 0),w x — —0.9 and three different dark 
energy densities: from bottom to top Q x — 0.7,0.74, and 0.8. Solid lines correspond to interacting models 
with c 2 = 10~ 3 ,0.01, and 0.1 (from bottom to top), w x = —0.9, and £l x = 0.74. Full diamonds represent 
measurements of the ISW amplitude obtained by cross-correlating different galaxy catalogs with WMAP first 
year data (see [4l| for details). For consistency with the data, the theoretical estimates are given in units of 
b. In the figure, the matter power spectra was normalized fixing the amplitude of the rms of matter density 
perturbations on a sphere of 8/i _1 Mpc to be a% = 0.8. As expected, if all models have the same as, the ISW 
effect increases with increasing fl x and increasing c 2 . The \ 2 P er degree of freedom is less than 1, except for 
c 2 = 0.01, fl x — 0.74 (Xdof — 1-5) an d c 2 = 0, Ojc = 0.8 (x^of — 3-6). An interaction parameter as high as 
c 2 = 0.1 is strongly ruled out by the observations, Xdof = 180. The low quality of the fit is dominated by the 



(T*M)(6) 



(10) 
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FIG. 4: Cross-correlation between CMB data and template, averaged on angular scales from 4 to 10°. Dashed lines 
corresponds to models with no interaction (c 2 = 0) and w x = —0.9 and, from bottom to top: Q x = 0.7,0.74, and 0.8. 
Solid lines correspond to interacting models (from bottom to top: c 2 = 10 _3 ,0.01, and 0.1), w x = —0.9, and Q x = 0.74. 
Full diamonds represent the measurements obtained from different catalogs and their la errors compiled in [4lj . 



measurement at z = 0.9 obtained in [5| and derived using NVSS and HEAO-1 samples. Due to the systematic 
uncertainties of the bias and selection function in both samples, our simple model may not be accurate enough 
to make a meaningfull prediction; but even if we disregard this data point, a value of c 2 = 0.1 can be ruled 
out at the 99.95% confidence level. The rejection level would be smaller if we considered the uncertainty in the 
shape of Um{z) and in the mean depth of galaxy catalogs. At any rate, this analysis shows that the ISW effect 
can effectively constrain the existence of an interaction at low redshift. It has more statistical power than the 
luminosity distance test based on SNIa observations [2lj ]. 

If WMAP first year data showed a good agreement between the best fit values of Sl x derived from the analysis of 
temperature anisotropy maps and those obtained by cross-correlating temperature maps with projected galaxy 
density templates, Cabre et al. (3fj| . |40T ] noticed some tension between the values derived from these techniques 
using the less noisy WMAP third year data. The ISW signature of the new CMB data had a preference for 
larger values of fl x . The combined analysis of different samples gave Q x = 0.77 — 0.86 at the 95% confidence 
level. The data was compatible with w x = — 1 (ACDM model), but did not have enough statistical power to 
break the degeneracy in the w x vs fl x plane. Those results were later confirmed in [391 ] where a best fit value 
of fl\ — 0.85 for a flat ACDM model was reported, albeit with larger uncertainty. An interaction parameter 
of c 2 ~ 0.01 could produce a ISW signal similar to a non-interacting model with higher fl x , bringing back the 
agreement between the measured fraction of dark energy with WMAP 3yr data and with the ISW effect. Should 
the Cabre et al. or Rassat et al. results [H, [H, H(| be confirmed by deeper catalogs like the new SDSS releases 
or PanSTARRS [44]], the ISW test will constitute an indirect evidence that the growth rate of matter density 
perturbations differs from that of the concordance ACDM model. Future surveys will be able to test different 
extensions of the standard model and prove/rule out the existence of an interaction. 



IV. CONCLUSIONS 



The ISW effect may provide an excellent tool to detect, in a model independent way, the coupling between dark 
energy and dark matter at low redshifts, i.e., z < 1. The interaction damps the growing mode of the Newtonian 
potential faster compared to models with no interaction and enhances the ISW effect; then, the ISW effect is 
sensitive to the growth rate of matter density perturbations, but model predictions need to be refined to model 
more realistically the selection function of galaxy catalogs. We have shown that the correlation of WMAP with 
templates constructed from different galaxy catalogs can rule out an interaction parameter of c 2 =0.1 at better 
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than the 99.95% confidence level. For this type of models, this upper limit is stronger than those obtained 
previously using the CMB radiation power spectrum [F9|, |2l|, but at a different redshift. The ISW effect shows 
more statistical power to discriminate among IQM models than the luminosity distance test obtained using 
SNIa. 

The recent results in (36l |40| suggest a discrepancy between the fraction of dark energy necessary to explain 
the amplitude of the ISW effect (measured from cross-correlating CMB data and galaxy catalogs) larger than 
allowed by the WMAP 3yr data. We have shown that the tension between both data sets could be explained 
by introducing a moderate coupling between dark energy and dark matter. This is in line with recent claims 
that the Layzer-Irvine equation [45| when applied to galaxy clusters reveals the existence of a DM/DE coupling 
. If the aforesaid discrepancy is confirmed by deeper and wider galaxy samples, it could turn out an indirect 
evidence for the existence of a DM/DE interaction. 
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